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Abstract. - Quantum spin transport is studied in an interacting quantum dot. It is found that a 
conductance "plateau" emerges in the non- linear charge conductance by a spin bias in the Kondo 
regime. The conductance plateau, as a complementary to the Kondo peak, originates from the 
strong electron correlation and exchange processes in the quantum dot, and can be regarded as 
one of the characteristics in quantum spin transport. 



Introduction. — The quantum charge transport has 
been extensively investigated in a wide class of correlated 
electron systems, which involves the interaction between 
a localized spin and free conduction electrons, such as in 
metals containing magnetic impurities [1]. Among them, 
the Kondo effect was of fundamental importance and was 
proposed to emerge in the device consisting of a single- 
electron quantum dot (QD) coupled to electrodes (i.e. 
leads) [2,3], and its conductance can be strongly enhanced 
in the Kondo regime due to the frequent occurrence of spin 
exchange processes. Rapid progresses in nano-technology 
have made it available to control the number of electrons 
in QD and the parameters (e.g. the energy level, coupling 
strength between QD and leads, and so on) of the QD [4]. 
The Kondo effect in the QD was observed experimen- 
tally [5,6], which was in excellent agreement with theoreti- 
cal predictions [2,3]. Most early works focus on the Kondo 
physics driven by the charge current, which was generated 
in two non-magnetic leads and under the bias voltage V 
(hereafter the bias voltage is named as charge bias). Re- 
cently the Kondo effect of the spin-polarized transport in 
QD with the ferromagnetic leads was investigated [7, 8] . 
It reveals novel properties of the Kondo effect in the spin 
aspect. Since electron has both charge and spin, it would 
be imcomplete for any physical picture of the Kondo ef- 
fect without full consideration of quantum spin transport. 



On the other hand, the spin current, in which electrons 
with spin up and down move in opposite directions, has 
attracted a lot of theoretical and experimental attentions. 
By using the magnetic tunneling injection, and electrical 
or optical injection technique, a pure spin current without 
accompanying the charge current can be generated and 
detected experimentally [9,10]. These experimental pro- 
gresses make it practical to investigate the properties of 
quantum spin transport in correlated systems. 

In the present paper, we introduce a spin bias instead 
of a charge bias to study the quantum transport in Kondo 
regime in an interacting QD system. The spin bias means 
the spin-dependent chemical potential in the two leads, 
and may produce a spin current flowing through the QD 
(see Fig. 2) [11]. Consider a common-used device, which 
consists of a QD coupled to two non-magnetic leads. It 
is found that a conductance plateau emerges in the non- 
linear charge conductance while a double peak appears in 
the non-linear spin conductance. The charge conductance 
plateau and spin conductance double peak only exist in 
the Kondo regime, and its width is equal to twice of the 
spin splitting of the QD's energy levels. It is believed that 
these phenomena can be regarded as one of the intrinsic 
characteristics in quantum spin transport. 

General formalism. 
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Model Hamiltonian and the spin bias. We start with 
the model Hamiltonian for the QD coupled to two leads [3] , 



Here 



H — HLead + Hp + Ht ■ 



H Lead = ]T e ak ci ka c aka 



(i) 

(2) 



describes the conduction electrons in the left (L) and right 
(R) leads a 6 (L,R) and c^ ak(y (c a ha) is the creation (an- 
nihilation) operator for electrons with momentum k and 
spin a E (t, 4) in the lead a. 



H D = ^e a n a + Un^n L (n a = d\d a ) 



(3) 



is the Hamiltonian for the interacting QD and, for sim- 
plicity, we only consider a single pair of levels with spin 
splitting e a = ed + erA/2. U represents the Coulomb in- 
teraction. 



(4) 



describes the tunneling between the QD and two leads 
with the tunneling coefficients t ak - The Kondo effect in 
this device has been extensively studied [2,3]. Most previ- 
ous works focused on the case of the charge bias V between 
two leads with the chemical potential fiLa = —^Ra = V/2, 
and part of the works were also concerned with the spin- 
polarized current driven Kondo effect [7, 8] . In the follow- 
ing, we apply the spin bias V, where 



Mil 



- m = V/2, 

-m = -v/2 



(5) 



(see Fig. 2) [11] to study the quantum electron transport. 
This spin bias can generate a spin current, in which the 
electrons with different spins move in opposite directions. 

Theory of Keldysh non- equilibrium Green's functions 
technique. In the theory of the Keldysh non-equilibrium 
Green's functions technique, the current from the lead a 
flowing into the QD for the spin channel a can be ex- 
pressed as [12], 

Iaa = liS ^ T ^Hf^(e) (G:(e)-G a a (e)) + G<(e)}, 

(6) 

where r a (e) = 2i:^ k t* ak t ak 5(e - e ak ) is the line-width 
function and / acr (e) = {exp [(e — /wV^b? 1 ] + 1} _1 is the 
Fermi-Dirac distribution of electrons in the lead a with 
spin a at temperature T. G^' a ' < (e) are the Fourier trans- 
formation of the retarded, advanced and lesser Green's 
function G r ^ a,< (t) in the QD, respectively, where G<(t) = 
i{dU0)d a (t)) andG^W =TiO(±t)({d a (t),dt(0)}). 

The Green's functions can be solved by the equation-of- 
motion technique, which was widely applied to study the 
electron transport in interacting QDs [3,13]. Although this 



method has some disadvantages in predicting the intensity 
of Kondo effect and the disappearance of the Kondo peak 
at ed + 2U = 0, it is relatively simple and straightforward, 
and has been proven to provide correct qualitative physics 
at low temperatures and in the large U limit (U — > oo). 
For the sake of simplicity, we concentrate on the transport 
properties in the strong interaction limit, i.e., U — > oo. In 
this case, the Fourier transform of the Green's functions 
G£ a (e) is obtained as, 



1 



S^(e)-S^(e) 



(7) 



where n a is the occupation number in the spin-er state 
and So' CT a (e) = 2tt l^afe | 2 /(e - e ak ± i0+) is the non- 



ah 



interacting tunneling self-energy. Sqct ( e ) can reduce to 
T,Q^(e) = ^(Tl + Tr)/2 while in the wide-band width 
limit, i.e. Tl and Tr being a constant independent of e. 
Si'^(e) in Eq.(7) due to the electron correlation and the 
tunneling has the form: 



cxk 



e - e ak -\- € & - € a ± i0+ 



(8) 



The occupation number n a in Eq.(7) should be solved self- 
consistently from the equation, n a = —i f(de/2n)G<(e). 

Next we come to solve the lesser Green's function G< (e) 
to calculate the current formula in Eq.(6) and the occupa- 
tion number n a self-consistently. For interacting systems, 
usually G<(e) cannot be solved exactly, and various ap- 
proximations were developed such as Ng ansatz [14], and 
the non-crossing approximation [3] . In the present calcula- 
tion, instead of G<(e) only J deG<(e) is actually used [15], 



J deG<{e) = -iJ delmGl(e) (f La {e) + f. 



R«(z))- (9) 



With the help of Eqs.(7), (8), and (9), we are ready to 
calculate the particle current I aa . Since there is no spin- 
flip process in the system, we have dn a /dt = —{Iro + Ilo) 
and I La = —Iru = la m a steady state. As a result, the 
charge current I c = e(/j + ij.) and the spin current I s = 
|(i| — I[) are conserved. In the following calculations, 
the line-width function r Q (e) is taken to be symmetric 
Tl(£) = T R (e) = T9(D - |e|), with a band width 2D = 
1000, and T = 1 as the unit of energy We define the 
charge differential conductance 



G C (V) = dI c /dV, 
and the spin differential conductance 
G S (V) = dI s /dV, 



(10) 



(11) 



with respect to the spin bias V. In the absence of the Zee- 
man splitting A = 0, if = — I[ because the Hamiltonian 
is invariant under the exchange of the indices (L, |) and 
(R, J,). In this case the charge current I c and its conduc- 
tance o~ c are equal to zero exactly. To obtain a non-zero 
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Fig. 1: (color online). The charge conductance dI c /dV in 
(a) and the spin conductance dI s /dV in (b) vs. the spin 
bias V at different temperatures T. The inset magnifies 
the curve of dI c /dV for k B T = 0.01. (T = 1, e d = -5, 
A = 0.4.) 

I c and <7 C the energy levels of the QD must be split, i.e. 
A ^ 0. For the spin current I s and its conductance a s they 
are always non-zero under a spin bias even if the energy 
levels are degenerate. 

It is worth pointing out that it is also possible to gen- 
erate a charge current as well as a spin current if the two 
tunneling coefficients are unequal, tt,jt fff.fr. even in the 
absence of the Zeeman splitting. The right-left asymme- 
try plus the strong Coulomb interaction may generate a 
charge current. The self-energy correction will also pro- 
duce the splitting of the energy level in the quantum dot, 
which looks like an effective "Zeeman" energy splitting, 
and varies with the strength of spin bias. However, the 
splitting is quite tiny. 

Results and discussions. — The main results in the 
present paper are summarized in Fig. 1(a) where a con- 
ductance plateau appears in the curve of the charge con- 
ductance dI c /dV and in Fig. 1(b) where a double peak 
appears in the curve of the spin conducatnee dI s /dV with 
respect to the spin bias V at different temperatures T. 
The double-peak occurs at V = ±A due to the energy 
splitting by a magnetic field. A similar result was ob- 
served from a solvable model. [16] It is also very simi- 
lar to the Kondo double-peak which is induced by the 
charge bias V c in the curve of the conventional charge 
conductance dI/dV c -V c , and has been observed experi- 
mentally [1,5,6]. The relationship between them will be 
discussed in more detail elsewhere. In the following, we 
focus mainly on the charge conductance induced by the 
spin bias. In Fig. 1(a), a plateau instead of the conven- 
tional Kondo-double-peak appears between —A < V < A, 



which is the key feature of spin-bias-induced conductance 
dI c /dV. This plateau is quite even, and the relative vari- 
ance \G C {V) - G c (0)|/G c (0) is within 0.1% in a quite 
large region of the spin bias voltage V (about 20fcsT^ 
at A = 0.4) [?]. At V — ±A, the conductance dI c /dV 
suddenly drops down within the scale of ksTx 1 • In par- 
ticular, this plateau has the same properties as the Kondo 
peak: it disappears completely if U — 0, or at higher 
temperatures T > Tk- The plateau survives only at low 
temperatures (T < Tk), and tends to saturate if T <C Tk 
(see Fig. 1(a)). 



a b 
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Fig. 2: (color online), (a) and (b) are the schematic di- 
agrams for the dominant exchange processes in \V\ < A 
and |V| > A, respectively. 

The conductance plateau is believed to originate from 
strong electron correlation and the exchange process as 
shown in Fig. 2. The strong Coulomb interaction only 
allows that the two levels in the QD can be occupied by 
at most a single electron. When \V\ < A as shown in Fig. 
2(a), there exist two exchange processes: (i) the spin- up 
electron in the QD tunnels into the right lead (the step 
1) and another spin-up electron in the left lead tunnels 
into the QD (the step 2) due to > £iRT> an d 00 the 
spin-down electron in the QD tunnels into the left lead 
(the step 1') and another spin-down electron in the right 
lead tunnels into the QD (the step 2') due to fXRi > [Xli- 
These two exchange processes may induce the charge cur- 
rent but with opposite flowing directions. Also because of 
Q| > e di at the positive Zeeman splitting A, the charge 
current from the exchange process (i) is larger than that 
from the process (ii), so that the conductance is positive 
and almost constant at \ V\ < A. On the other hand, while 
|V| > A, the spin exchange processes as shown in Fig. 
2(b) also occur. However, these exchange processes are 
prohibited when \V\ < A, because the energy is required 
to be conserved after the two virtual tunneling steps 1 
and 2 , and 1' and 2' in Fig. 2(b). In particular, these 

1 Although the Kondo temperature Tk = D exp[— it^d/i^L + 
T R )\ »s 0.2 for e d = -5T, T K is estimated to be about 0.04 
from the half-width of the Kondo peak in Fig. lb. Usually the 
cquation-of-motion method underestimates the electron correlation 
and gives a lower value of Tk- For the discussion in the text, we 
take k B T K = 0.04 for e d = -5T. 
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exchange processes can frequently occur at low temper- 
atures, so that the exchange processes in Fig. 2(a) are 
restrained. As a result the conductance dI c /dV sharply 
drops while |V| passing through A. As for the weak in- 
teraction limit, U = 0, the QD can be occupied by two 
electrons simultaneously, and the exchange processes in 
Fig. 2(b) do not occur because of the Pauli exclusion 
principle, and the conductance pleatau disappears. The 
co-tunneling processes in Fig. 2(b) also explain the occur- 
rence of the double-peak in the spin conductance dI s /dV. 
These processes in step 1 and 2 , and 1' and 2' exchange 
the spins in the two spin channels in the leads. At each 
process of step 1 and 2 , and 1' and 2' it will produce spin 
current instead of charge current. Thus a resonance of 
spin differential condctance emerges at the critical points 
V = ±A and the double-peak arises correspondingly. 




| A=0.0 
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Fig. 3: (color online). The charge conductance dI c /dV 
vs. the spin bias V for different A. (fcsT = 0.01 and 
e d = -5.) 

To explore the physical properties of the plateau in 
charge conductance, we come to study the variance of 
the plateau with A and ed- Fig. 3 shows the conduc- 
tance dI c /dV versus the spin bias V for different A at 
T < T K and e d = -5. When A = 0, dI c /dV = for 
any V and no plateau appears. With increasing A from 
zero, dI c /dV\ v=0 rises very quickly. The plateau begins 
to form when A is about a few ksTK- Its width is 2 A and 
its height depends on A weakly. This plateau remains 
robust even when A reaches about 20fcsTR-. At a larger 
A, a slightly downward bend emerges in the center of the 
plateau. 

Fig. 4 shows the e^-dependence of the conductance 
dI c /dV versus the spin bias V at T < Tr. The conduc- 
tance plateau always exists and is robust when ed < —5, 
i.e. in the Kondo regime where the electron is localized in 
the QD. With increasing ed, the electron is partially delo- 
calized leading to the plateau distorted. When the system 
is in the mixed- valence regime, i.e., —5 < < 0, the 
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Fig. 4: (color online). The charge conductance dI c /dV vs. 
the spin bias V for different e d . The inset plots the curve 
dI c /dV-ed at V = 0. The arrows indicate the values of ed 
in the main figure. (A = 0.4 and k B T = 0.01.) 

plateau disappears completely, but the two sharp jumps 
still survive at V — ±A, and a deep valley emerges at 
V = 0. In fact, the linear charge conductance G c (0) is 
intensively dependent on in the mixed-valence regime. 
The inset of Fig. 4 plots G c (0) versus ed, in which a deep 
dip exhibits at about e d = — 2 and G c (0) can even change 
its sign. At last, when > 0, the QD is in an empty 
state, and the dI c /dV-V curve becomes relatively smooth 
due to the nonexistence of the exchange process at ed > 0. 
This behavior suggests that the conductance plateau in 
the Kondo regime is closely related to the localization of 
electron in the Kondo regime, and the physics behind the 
plateau deserves further discussion. 

Finally we come to discuss the combined effect induced 
by the spin and charge bias. Consider the spin bias V s = V 
and the charge bias V c — Vq, i.e., — ±V/2 + Vo and 

/| = =p V/2. In this case both spin and charge current 
through the QD are non-zero. Fig. 5 plots the dI c /dV-V 
curve for different charge bias voltages Vo- The results 
exhibit very well the combination of the plateau and the 
peak. The peaks are of Kondo type induced by charge 
bias. For a small Vq (see Fig. 5(a)), in addition to the 
plateau, two extra Kondo peaks and dips begin to emerge 
in the dI c /dV-V curve at the positions V — ±A ± Vo- 
The Kondo dips originate from the fact that the voltage 
difference — jiR^ decreases linearly with increasing of 
the charge bias Vo, and the conventional Kondo peaks from 
the spin-down electron are transformed into the dips. On 
the other hand, for a large Vo (see Fig. 5(b)), the Kondo 
peaks and dips become dominant. 

Summary. — In summary, a conductance plateau 
emerges in the charge differential conductance dI c /dV-V 
curve in quantum spin transport, instead of the Kondo- 
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Fig. 5: (color online), (a) and (b) show the curve dI c /dV- 
V while both spin and charge bias coexist, (e<j = —5, 
A = 0.4, and k B T = 0.01.) 



double peak in quantum charge transport. The conduc- 
tance plateau and Kondo peak are complementary to each 
other, and reflect the spin and charge aspects of the ex- 
change processes in quantum electron transport. 
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